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The magnetization of the prototypical molecular magnet Mni2-acetate 
exhibits a series of sharp steps at low temperatures due to quantum 
tunneling at specific resonant values of magnetic field applied along the 
easy c-axis. An abrupt reversal of the magnetic moment of such a crystal 
can also occur as an avalanche, where the spin reversal proceeds along a 
"deflagration" front that travels through the sample at subsonic speed. 
In this article we review experimental results that have been obtained 
for the ignition temperature and the speed of propagation of magnetic 
avalanches in molecular nanomagnets. Fits of the data with the theory of 
magnetic deflagration yield overall qualitative agreement. However, nu- 
merical discrepancies indicate that our understanding of these avalanches 
is incomplete. 

1.1. Background 

First synthesized by Lis in 1980f3 Mni 2 Oi2(CH3COO)i 6 (H 2 0)4 (referred 
to hereafter as Mni2-ac) is a particularly simple, prototypical molecular 
magnet. Shown in Fig- ILlf a). the magnetic core of Mni2-ac has four Mn 4+ 
(S = 3/2) ions in a central tetrahedron surrounded by eight Mn 3+ (S = 2) 
ions. The ions are coupled by superexchange through oxygen bridges with 
the net result that the four inner and eight outer ions point in opposite 
directions, yielding a total spin S — 10 P The magnetic core is surrounded 
by acetate ligands, which serve to isolate each core from its neighbors in a 
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body-centered tetragonal lattice. A crystalline sample contains ~ 10 17 or 
more (nearly) identical, weakly interacting single molecule nanomagnets in 
(nearly) identical crystalline environments. 

While there are very weak exchange interactions between molecules, the 
exchange between ions within the magnetic core is very strong, resulting 
in a rigid spin 10 unit that has no internal spin degrees of freedom at low 
temperatures. To lowest order, the spin Hamiltonian is given by: 



The first term denotes the anisotropy barrier, the second is the Zeeman 
energy that splits the spin-up and spin-down states in a magnetic field, and 
the last term, W ', contains all symmetry- breaking operators that do not 
commute with S z . For Mni2-ac, D = 0.548K, g z = 1.94, and /is is the 
Bohr magneton. 



Fig. 1.1. (a) Chemical structure of the core of the Mni2 molecule. The four inner 
spin-down Mn 3+ ions each have spin S = 3/2; the eight outer spin-up Mn 4+ ions each 
have spin 5 = 2, yielding a net spin S = 10 for the magnetic cluster; the small open 
circles are O bridges and arrows denote spin. Acetate ligands and water molecules have 
been removed for clarity, (b) Left: Double-well potential in the absence of magnetic field 
showing spin-up and spin-down levels separated by the anisotropy barrier. Different spin 
projection states \m > are indicated. The arrows denote quantum tunneling. Right: 
Double-well potential for the N=2 step in a magnetic field applied along the easy axis. 

As illustrated by Figure ITTlT b). the spin's energy can be modeled as a 
double- well potential, where one well corresponds to the spin pointing "up" 
and the other to the spin pointing "down". A strong uniaxial anisotropy 
barrier of the order of 66 K yields doubly degenerate ground states in 
zero field. The spin has a set of energy levels corresponding to different 
projections, m = 10, 9, . . . , — 9, — 10, of the total spin along the easy (c- 
axis) of the crystalPS 



U = -DSl - g zf i B H z S z + ...+H'. 



(1.1) 
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1.2. Spin Reversal by Quantum Tunneling 

Slow relaxation below the blocking temperature, Tb ~ 3 K, gives rise to hys- 
teresis loops that display steps^as a function of magnetic field, H z , swept 
along the easy c-axis of the crystalP^ Figure fL2l a) shows the magnetiza- 
tion M as a function of magnetic field HqH z \ the derivative, dM /dH, which 
reflects the magnetic relaxation rate, is plotted as a function of \x§H z in 
Fig. ll.2f b). These steps, characteristic of molecular magnets, can be under- 
stood with reference to the double well potential of Fig. Il.lf b): a magnetic 
field H z introduces a Zeeman splitting that tilts the potential wells and 
causes energy levels in the right (left) well to move down (up). Levels in 
opposite wells align at particular values of magnetic field (dashed lines in 
Fig. Il.lf b)). allowing the spin to reverse by tunneling. Full (saturation) 
magnetization is thereby reached in a stepwise fashion, with the detailed 
form of the steps depending on sweep-rate and temperature. 

1.3. Spin Reversal by Avalanches 

As first reported by Paulsen and Parkf^l Mni2-ac crystals sometimes ex- 
hibit a sudden, complete reversal of magnetic moment during a field-swept 
measurement. This phenomenon, also observed in other molecular magnets, 
was attributed to a thermal runaway (avalanche) in which the relaxation 
of magnetization toward the direction of the field results in the release of 
heat that further accelerates the magnetic relaxation. Direct measurements 
of the heat emitted have confirmed the thermal nature of the avalanches. 
In addition to releasing thermal energy, molecular crystals emit bursts of 
radiation during magnetic avalanches EHED Once considered events to be 
avoided, as they interfere with a detailed study of the stepwise process of 
magnetization, magnetic avalanches became the focus of attention and re- 
newed interest stimulated by the theoretical suggestion that the radiation 
emitted during an avalanche is in the form of coherent (Dicke) superra- 
dianceP^I Although the issue of coherence of the radiation has yet to be 
resolved, recent studies have clarified the nature of the avalanche process 
itself. 

1.3.1. Magnetic Deflagration 

From time-resolved measurements of the local magnetization using an ar- 
ray of micron-sized Hall sensors placed on the surface of Mni2-ac crystals, 
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Fig. 1.2. (a) Magnetization of a Mni2-ac crystal normalized by its saturation value 
as a function of magnetic field applied along the uniaxial c-axis direction at different 
temperatures below the blocking temperature; the magnetic field was swept at 10 mT/s. 
(b) The derivative, dM/dH of the data in part (a) as a function of magnetic field. 




Fig. 1.3. (a) The local magnetization measured as a function of time by an array of 
micron-sized Hall sensors placed along the surface of the sample. The inset illustrates 
the "bunching" of magnetic field lines as the deflagration front travels past a given Hall 
sensor. (b)The sensor position as a function of the time at which the sensor registered the 
peak; the inset shows the placement of the Hall sensors on the crystal. The propagation 
speed for this avalanche is 2.2 m/s, approximately three orders of magnitude below the 
speed of sound. 



Suzuki et aZP^ discovered that a magnetic avalanche propagates through 
the crystal at subsonic speed in the form of a thin interface between regions 
of opposite magnetization. Figure [T.3.1f a) shows traces recorded during an 
avalanche by sensors placed in sequential positions near the center of a 
Mni2-ac sample. The inset is a schematic that illustrates the bunching of 
field lines at the propagating front that gives rise to the observed peaks. 
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Figure 11.3.1( b) is a plot of the sensor position versus the time of arrival of 
the peak. From these measurements one deduces that the front separating 
up- and down-spins travels with a constant (field-dependent) speed on the 
order of 5 m/s, two to three orders of magnitude slower than the speed of 
sound. 

From a thermodynamic point of view, a crystal of Mni2 molecules placed 
in a magnetic field opposite to the magnetic moment is equivalent to a 
metastablc (flammable) chemical substance. A well-known mechanism for 
the release of energy by a metastable chemical substance is combustion 
or slow burning, technically referred to as deflagration!^ It occurs as a 
flame front of finite width propagates at a constant speed small compared 
to the speed of sound. For "magnetic deflagration" in Mni2-ac, the role of 
the chemical energy stored in a molecule is played by the difference in the 
Zeeman energy, AE = 2gfiBHS, for states of the Mni2-ac molecule that 
correspond to S parallel and antiparallel to H. 

1.3.2. Avalanche Ignition 

Although the probability of a spontaneous avalanche has been shown to 
be higher at resonant magnetic fields than off-resonancep' avalanche igni- 
tion is unpredictable and uncontrolled when an external magnetic field is 
swept back and forth, the experimental protocol generally used to study the 
steps in the hysteresis loops. Avalanche ignition under these conditions is a 
stochastic process that depends on factors such as the sweep rate, the tem- 
perature, and the quality of the crystal. Controlled ignition of avalanches 
has now been achieved using surface acoustic waves (which serve to heat 
the sample) and by using a heater^ , as described below and in the next 
section. 

McHugh et al}® employed a resistive wire element as a simple electric 
heater to trigger avalanches in a manner similar to the work of Paulsen and 
ParkP^l In these experiments, an external magnetic field is ramped to and 
held at a fixed value. The heater is then turned on to slowly heat the sam- 
ple until an avalanche is triggered at a temperature measured by a small 
thermometer. Avalanches launched by this method occur at well-defined, 
reproducible ignition temperatures. Fig. 11.4( a) shows a typical tempera- 
ture profile: starting at the base temperature of 300 mK, the temperature 
gradually rises until an abrupt sharp increase in the temperature signals 
the ignition of an avalanche. For this avalanche triggered at [ioH z = 1.85 
T, the ignition temperature is about 0.6 K. 
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Fig. 1.4. (a) Temperature recorded by a thermometer in contact with a Mni2 crystal 
during the triggering of an avalanche. The heater is turned on at ~ 0.012 s, the tem- 
perature then increases slowly until an abrupt rise in temperature at 0.03 s signals the 
ignition of an avalanche. The inset shows data taken near ignition with higher resolu- 
tion. The noise at low temperatures derives from digitizing the analog output of the 
thermometer, which depends weakly on temperature below 0.4 K. (b) Temperature pro- 
files for avalanches of major and minor species triggered at low fields in a Mni2 crystal. 
The two types of avalanches are triggered separately below a sample-dependent magnetic 
field, while at higher fields ignition of the minor species triggers the ignition of the major 
species. 

Single crystals of M1112— ac are known to contain two types of molecules. 
In addition to the primary or "major" species described earlier, as-grown 
crystals contain a second "minor" species at a level of w 5 percent with 
lower (rhombohedral) symmetry! 20 * 21 ! These faster-relaxing molecules can 
be modeled by the same effective spin Hamiltonian, Eq. 11.11 with a lower 
anisotropy barrier of 0.49 K. Avalanches of the each species can be stud- 
ied in the absence of the other through an appropriate magnetic protocol 
described in Ref. [22] Interestingly, avalanches are separately triggered 
by the two species in low magnetic field. As shown in Fig. 11.4( b). at low 
fields the minor species relaxes prior to and independently of the major 
species, while above ~ 0.7 T, the major and minor species ignite together 
and propagate as a single front. It is analogous to grass and trees that can 
sustain separate burn fronts that abruptly merge into a single front when 
the grass becomes sufficiently hot to ignite the trees. 

Despite the turbulent conditions that one might expect for deflagration 
(analogous to chemical combustion), quantum mechanical tunneling clearly 
plays a role, as demonstrated in Fig. 11.51 where the temperature above 
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Fig. 1.5. Temperature required to ignite avalanches plotted as a function of magnetic 
field. The vertical lines denote the magnetic fields where sharp minima occur in the 
ignition temperature corresponding to tunneling near the top of the anisotropy barrier. 
The overall decrease in ignition temperature is due to the reduction of the anisotropy 
barrier as the field is increased. 

which ignition occurs is plotted as a function of a preset, constant magnetic 
fieldP^ The temperature required to ignite avalanches exhibits an overall 
decrease with applied magnetic field, reflecting the fact that larger fields 
reduce the barrier (see the double- well potential in Fig. Il.lf b)). The role 
of quantum mechanics is clearly evidenced by the minima observed for the 
ignition temperature at the resonant magnetic fields due to tunneling when 
levels cross on opposite sides of the anisotropy barrier. 

In the ignition studies described above, the barrier against spin reversal 
was lowered by applying an external magnetic field, H z , along the uniaxial 
c-direction, which serves to unbalance the potential wells and lower the bar- 
rier to spin reversal. A transverse field, H x , also reduces the anisotropy bar- 
rier by introducing a symmetry-breaking term, {g^BH x S x ), to the Hamilto- 
nian, Eq. 11.11 thereby promoting tunneling. Macia et alW^ investigated the 
threshold for avalanche ignition in Mni2 -0jC clS cl function of the magnitude 
and direction of a magnetic field applied at various angles with respect to 
the anisotropy axis and as a function of temperature. As the external field 
is increased at a constant rate from negative saturation to positive values, 
both H z and H x increase, tracing a trajectory in the (H Z ,H X ) parameter 
space; (examples of sweeps starting from zero are shown by the arrows in 
Fig. 11.61 As shown in Fig. 11.61 fa), an avalanche was recorded for each pair 
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Fig. 1.6. (a) Angle dependence of metastability measured through the occurrence of 
avalanches. Squares (triangles) denote parameter values where deflagration occurs for 
initial temperature 2.2 K (1.8 K). (b) Theoretical calculation for the area of stability 
against ignition of avalanches (solid curve) and against slow relaxation (dashed curve). 
Circles denote points where avalanches are predicted to occur at a given angle within 
the first quadrant. The angle 8 C denotes the crossing point between areas of slow re- 
laxation and avalanche stability. These results were obtained with as a parameter 
varying from 6.8 K for H = 4600 Oe to 10.9 K for H = 9200 Oe. From Macia et aZ.GH 



(H X ,H Z ) denoted by a square (for T = 2.2 K) or a triangle (for T = 1.8) 
K. We postpone a detailed explanation of these results and a comparison 
with theory to a later section. 



1.3.3. Avalanche Speed 

Following the initial studies of Suzuki et aZP^ in which avalanches were trig- 



gered stochastically in swept magnetic field, Hernandez-Minguez et 
carried out a systematic investigation of avalanche speeds as a function of a 
preset, constant magnetic field, HoH z , for avalanches triggered controllably 
using surface acoustic waves. From SQUID-based measurements of the to- 
tal magnetization of a crystal of known dimensions, and the realization 
that the avalanche propagates as an interface between regions of oppo- 
site magnetization,^ Hernandez-Minguez et a?! 18 * 23 ! deduced the speed of 
propagation as a function of magnetic field shown in Fig. 11.71 The speed 
of the avalanches is enhanced at the resonant fields where tunneling oc- 
curs, confirming the important role of quantum mechanics, and prompting 
the authors to name the phenomenon "quantum magnetic deflagration." 
Similar results have been obtained from local, time-resolved magnetization 
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Fig. 1.7. The speed of propagation of the magnetic avalanche deflagration front is plot- 
ted as a function of the (fixed) field at which the avalanche is triggered; note the enhance- 
ment of propagation velocity at magnetic fields corresponding to quantum tunneling 
(denoted by vertical dotted lines). From Hernandez-Mi'nguez et aZ.E^l 

measurements using micron-sized Hall sensors. 

McHugh et alpS have reported a detailed, systematic investigation 
of the speed of magnetic avalanches for various experimental conditions. 
The speed of propagation of an avalanche is described approximately^ by 
the expression, v ~ (K/ro) 1 ^ 2 exp[— U(H)/2kBTf], where U is the barrier 
against spin reversal, Tf is the flame temperature at or near the propagat- 
ing front where energy is released by the reversing spins, k is the thermal 
diffusivity, and tq is an attempt time. In these studies, avalanches were 
controllably triggered: (A) in various external fields with fixed (maximum) 
initial magnetization, so that both U and Tf vary; (B) in fixed external field 
with different initial magnetization so that the avalanches differ primarily 
through the amount of energy released, and thus the flame temperature 
Tf; and (C) where external magnetic fields and initial magnetization are 
varied and adjusted to hold the energy released, and thus Tf, constant. 
These parameters will be discussed in more detail in the next section. We 
note at this point that the energy barrier U and the flame temperature Tf 
appear only as the ratio U /Tf in the above expression for the velocity. It 
is therefore convenient to plot the speed of the avalanche as a function of 
U/Tf as is done in Fig. Ota)P 

To end this section, we note that Villuendas et alP*l have recently re- 
ported a novel approach to studying avalanche dynamics using the magneto- 
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Fig. 1.8. (a) Avalanche speeds for a single crystal with various initial magnetic prepa- 
rations. A denotes avalanches with A.M/2M S = 1; B\ and B2 denote data taken at 
(J-oHz = 2.2 T and 2.5 T, respectively; C\ and C2 denote avalanches with estimated 
flame temperatures Tf m 10 K and 12 K, respectively, (b) Avalanche speeds for different 
crystal with AM/2M S = 1. The fit requires an unphysical temperature dependence for 
the thermal diffusivity, k oc T 3 ' 5 . 



optical Kerr effect to image the magnetization during an avalanche. Al- 
though the resolution was quite limited, with further improvements this 
technique could provide a valuable avenue for studying magnetic deflagra- 
tion. 



1.4. Comparison with the Theory of Magnetic Deflagration 

There are two essential ingredients for magnetic deflagration: the 
metastable spins release energy as they relax to the ground state, and this 
energy diffuses as heat through the crystal and thermally stimulates the 
reversal of neighboring metastable spins. Garanin and Chudnovskyl^H de- 
veloped a comprehensive theory of magnetic deflagration encapsulated in 
the following equations!^ 

c !g_v.*vr = -( E )f (1.2) 

£-1* (1.3) 

Equation 11.21 describes the flow of heat through the solid with the relaxing 
spins, ^jr, as a source of heat; the thermal conductivity, k, is related to the 
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heat capacity, C, through the thermal diffusivity, k, as k — kC . Equation 
11.31 describes the relaxation of the metastable spin density, n, with a ther- 
mal relaxation rate T given by an Arrhenius law, r = Tq exp[— U /fesT]; (E) 
is the average amount of heat released per molecule when its spin relaxes to 
the stable state. For a single relaxing molecule, the energy change from the 
metastable to stable state is simply the Zeeman energy AE — 2gfisSB z . 
However, since not all spins necessarily relax during an avalanche, the av- 
erage energy released per molecule must be introduced: 



where M s is the saturation magnetization and AM — \M Z — M s \ is the 
change from initial to final magnetization. 

Since both k and C are functions of temperature, Eqs. 11.21 and 11.31 are 
coupled, nonlinear differential equations, which make them very difficult 
to solve in general. The samples used in the experiments generally have 
large aspect ratios, so that we can use simple approximations for quasi-one- 
dimensional avalanche propagation 26 to compare theory with experimental 



We begin by comparing the theory with experimental results for the 
ignition temperatures and the stability criteria. Heuristically, a deflagration 
front can develop when the rate at which energy is released by the relaxing 
metastable spins exceeds the rate of energy loss through the boundaries of 
the crystal. This condition can be expressed in terms of a critical relaxation 



where T is the initial temperature, T c — T exp[— [//fesTb], and I 2 is the 
characteristic cross section of the crystal. The curves shown in Fig. 11.61 (b) 
are the result of a calculation using Eq. 11.51 Two areas are defined in the 
(H Z ,H X ) parameter space where the spins are expected to be metastable 
against relaxation: the solid line denotes the region of metastability against 
relaxation by triggering avalanches while the dashed curve delineates the 
region of metastability against slow, stepwise relaxation. If the experimen- 
tal trajectory, denoted by the arrows, crosses the grey solid line first, an 
avalanche will ignite. If the dashed line is crossed first, the metastable spins 
will relax slowly without triggering an avalanche. This defines a critical an- 
gle 9 C , above which an avalanche cannot occur. 




(1.4) 



results. 



ratepH 



8k(T Q )k B T 2 
U(E)l 2 



(1.5) 
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Macia et al. measured the ignition threshold by applying an increasing 
external field at an angle with respect to the crystal. The relaxation rate 
increases as the field grows until r c is reached and deflagration ignites, as 
shown in Fig. 11.61 (a). For sufficiently large values of H x , they found that 
the slow relaxation of the metastable spins occurs before deflagration can 
ignite. This defines a line in parameter space separating regions where one 
or the other mode of relaxation occurs, as shown in Fig. I1.6f b). The theory 
predicts that the transverse field should result in a significant decrease in 
the magnetization metastability at the resonant fields of H z . The data 
recorded in Fig. I1.6f a) confirm this and are consistent with the ignition 
temperatures of Fig. 11.51 In addition, ignition thresholds were measured 
at two different temperatures. The area of stability is clearly reduced by 
the increased initial temperature, as expected. 

To summarize, the theory of magnetic deflagration is in excellent agree- 
ment with the experiments of McHugh et al., where the critical relaxation 
rate was reached by varying To with a heater, and the experiments of Macia 
et at, where the ignition threshold was reached by controlling the barrier 
U using H x and H z . 

We now proceed to compare the measured avalanche speeds with the 
theory of magnetic deflagration. An approximate expression for the speed 
of the deflagration front is given by^fil 



Eq. I1.6l requires that the relaxation rate of the spins at the highest temper- 
ature, T(U, Tf), is significantly slower than the rate at which heat traverses 
the interface width, k/5 2 , where S is the magnetic interface width; this 
condition is expected to hold for all speeds considered here. 

The barrier U{B Z = ^oH z ) is calculated from the effective spin Hamil- 
tonian (Eq. II. 1[) . The temperature of the front, Tf, can be estimated from 
the measured heat capacity at zero magnetic field the calculated contri- 
bution of the spins in magnetic field, and the average energy released per 
molecule: 



J T 

Typical values calculated for Tf range from 7 K to 16 K. 

With reference to the theoretical expression for the avalanche speed, 
Eq. 11.61 if one assumes the thermal diffusivity k is approximately indepen- 




(1.6) 




(1.7) 
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dent of temperature, or that its temperature dependence is unimportant 
compared to that of other parameters in the problem, then all measured 
avalanche velocities should collapse onto a single curve when plotted as a 
function of [U(H)/Tf]. Although an approximate collapse is indeed ob- 
tained, as shown in Fig. 11.81 (a), there are clear and systematic deviations 
depending on whether (A): there is full (maximum) magnetization rever- 
sal, AM/2M S = 1 (both U and T f vary); (Bi,B 2 ): the amount of "fuel" 
AM/2M is varied for a fixed magnetic field (thus U is held constant); 
(Ci,C2): avalanches are triggered such that the product AM x H remains 
the same so that the energy released and the flame temperature Tf are held 
constant. That these different experimental protocols introduce systematic 
variations, albeit small, suggests that the theory is incomplete. 

An attempt to fit to the theory by allowing the thermal diffusivity to 
vary with temperature as a power law is shown in Fig. 11.81 (b) for avalanches 
of type (A) that involve full magnetization reversal. Note that the enhance- 
ments of the velocity at certain values of magnetic field are associated 
with the tunneling resonances, in agreement with Fig |1.7l The fit with Eq. 
11.61 is obtained for a thermal diffusivity that varies with temperature as 
K ~ T 3 - 5 . This is distinctly unphysical, as the thermal diffusivity is gen- 
erally a strongly decreasing function of temperature^ for these materials. 
A similar analysis was performed^ on avalanches of the minor species that 
yields an even steeper increase of ft with temperature. Regrettably, exper- 
imental measurement of the thermal diffusivity of Mni2 are not available. 

We conclude that, although the theory of deflagration agrees well with 
the measured conditions for the ignition of avalanches and provides a de- 
scription of the avalanche velocity that is qualitatively correct, there are de- 
tailed discrepancies that suggest that additional factors need to be included 
in the theory to obtain good quantitative agreement with experiment. 

1.5. Outlook 

More experimental work is clearly needed. Measurements of the thermal 
diffusivity would provide an important constraint on the theory, as would a 
reliable determination of the (local) temperature of the deflagration front. 
Investigations of the influence of sample shape, size and quality would also 
be useful. Spatial control of the avalanche ignition points, possibly by 
optical means, could provide important information. Studies of the shape 
of the deflagration front, and its character (turbulent or laminar) would be 
particularly interesting. 
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The possibility of observing detonation is intriguing. Deflagration is but 
one type of combustion process. Another, more violent type, is detonation, 
where heat spreads from the reaction front as a shock wave rather than by 
diffusion. It is natural to ask whether crystals of molecular magnets can 
support the magnetic analog of chemical detonation. Decelle et alP^ have 
reported intriguing results hinting at this possibility using high external 
field sweep rates (4 kT/s). The interpretation of these experiments is not 
entirely clear, and much work remains to be done. 

We close by noting that, to the degree that magnetic deflagration re- 
sembles chemical deflagration, the magnetic manifestation of this process 
offers some major advantages for the systematic study of chemical combus- 
tion. The magnetic analog is non-destructive, reversible and continuously 
tunable using an external field. Unlike the chemical process, it is a partic- 
ularly interesting realization of deflagration in which quantum mechanical 
tunneling plays an important role. 
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